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1. Introduction 

The excitations of the QCD vacuum with different quantum numbers as well as their 
life times and decay modes are of fundamental interest in the physics of the strong 
interaction. The masses of the pseudo-scalar mesons have been found to be essentially 
due to a spontaneous breaking of the chiral SU(3)r x SU(3)l symmetry or the U(1)a 
anomaly (in case of the rf). The vector mesons p, u, 0, K*, J/^S etc., which are the 
dipole modes of the vacuum, have found increasing attention during the last two decades. 
Especially their decay to dileptons is presently investigated in elementary and complex 
(nucleus-nucleus) collisions in different laboratories all over the world (cf. the reviews 
IH|, |2|, |3| and Refs. therein). On the other hand, the scalar sector of vacuum excitations 
is not well known experimentally and theoretically, so far. 

The structure of the lightest scalar mesons ao(980) and /o(980) is still under 
discussion (see e.g. [|], ^ ||, [7|, [8], B anc ^ references therein). Different authors 
interpreted them as unitarized qq states or as four-quark cryptoexotic states or as KK 
molecules or even as vacuum scalars (Gribov's minions). Although it has been possible 
to describe them as ordinary gg-states (see Refs. |TT], [12], |i~3fl), other options cannot 
be ruled out up to now. Another problem is the possible strong mixing between the 
uncharged a (980) and the /o(980) due to a common coupling to KK intermediate states 



15| , |16| , [H]],[]n], [L9], p(| . This effect can influence the structure of the uncharged 
component of the ao(980) and implies that it is important to perform a comparative 
study of a(] and a$ (or a^). There is no doubt that new data on a° and a$ /oq production 
in iiN and NN reactions are quite important to shed new light on the Oq structure and 
the dynamics of its production. 

In our recent paper |2l] we have considered ao production in the reaction ttN — > a^N 
near the threshold and at GeV energies. An effective Lagrangian approach as well as 
the Regge pole model were applied to investigate different contributions to the cross- 
section of the reaction 7rN —>■ a N. Here we employ the latter results for an analysis 
of ao production in NN collisions. Our study is particularly relevant to the current 
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Our paper is organized as follows: In Sect. 2 we discuss an effective Lagrangian 
approach with one-pion exchange while the Reggeon exchange model is considered in 
Sect. 3. Sect. 4 is devoted to the calculations of the cross section for the reaction 
NN — > NNoq. In Sect. 5 we analyze the contribution of the ao resonance to the cross 
sections and invariant KK mass distributions for the reactions pp — > ppK + K~ and 
pp —>■ pnK + K°. Our conclusions are presented in Sect. 6. 

2. An effective Lagrangian approach with one-pion exchange 

We consider a{], aj", a^ production in the reactions j = pp — > ppa^ pp -^ pna^, 
pn — > ppaQ and pn — > pna^ using the effective Lagrangian approach with one-pion 
exchange (OPE). For the elementary nN —* Noq transition amplitude we take into 
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account different mechanisms a corresponding to t-channel diagrams with r/(550)- and 
/i(1285)-meson exchanges (a = t(n), £(/i)) as well as s- and u-channel graphs with an 
intermediate nucleon (a = s(N), u(N)) (cf. Ref. pifl). The corresponding diagrams 
are shown in Fig. [T| The invariant amplitude of the NN — > NNoq reaction then is the 
sum of the four basic terms (diagrams in Fig. [I]) with permutations of nucleons in the 
initial and final states 

M] (a) [ab; cd] = q {a) [ab- cd] M* a [ab\ cd] + % {a) [ab; dc] Ml[ab; dc] (1) 

+ C j{a) [ba- dc] Ml[ba- dc] + £] {a) [ba; cd] M^ba; cd], 

where the coefficients £T/ a ) are given in Table |l|. The amplitude for the t-channel 
exchange with ??(550)- and /i(1285)-mesons are given by 

Mt {v) [ab; cd] = g aoV nF aovn ((p a - p c ) 2 } (p d - p b f) g n NNF v ((p a - p c ) 2 ) 

x 7 ^2 2 ^(PchMPa) x Tl(p b ;p d ), (2) 

M\ h) [ab; cd] = - ga f^F aofl7T ((p a - p c ) 2 , (p d - p b ) 2 ) g flNN F fl ((p a - Pc. 

1 / . o / \\ ( (Pa-Pc)n(Pa-Pc)i 



(Pa -Pc) 2 - m 2 fl 



(Pa-Pc + 2 [Pb-Pd))n [gyv 



9 \fa fc i " \fo fa a I 3uv 9 

\ m 



h 



x u(p c )j 5 j v u(p a ) x Il(p b ;p d ), (3) 

with 

n(Pb]Pd) = -^ i^ ((P6 -Pd) 2 ) (Pb ~ Pd)f}u(pdh5l/pu(pb) 7 r^ 2- ( 4 ) 

The amplitudes for the s- and u-channels (lower part of Fig. |I]) are given as 

\a-k r u J] ni \ f* NN t? (t \2\ F N ((p a + p b - p d ) ) 

M s(N) [ab;cd]=U(p b ;p d ) —F„{{p d -p b ) ) g aoNN ^ — _ -^ _ -^ 

x (Pd - P&) M «(Pc) [(Pa + Pb - Pd)<57<5 + m N ]-f 5 j^u(p a ) , (5) 

A^vr r 7 ,1 t-t/ \ f-wNN j? (1 \2\ F N ((p c + p d - p b ) ) 

M u{N) [ab-cd] = U(p b]Pd ) —F„ {( Pd - Pb ) ) g aoNN - — — _ -^ 

x (Pd - P&)m w(Pc)757m[(Pc + Pd - P&)<57<5 + mjv]u(Pa)- (6) 

Here p a , p?, and p c , p d are the four momenta of the initial and final nucleons, respectively. 
The effective Lagrangians involving a$ and f\ mesons were taken in the following 
forms: 

^aorjn 9aorfK i]\X) 71 [X ) Qq\X J , 

£a /i7T = £a /ivT ef(x) <9 A 7r(x) O (x), (7) 

£a JV7V = g ao NN *jv(x) a (x) ^n(x), 

£fiNN = gfiNN e{'(x) ^ N (x) 7 A *iv(x). 

We mostly employ coupling constants and form factors from the Bonn-Julich potentials 
(see e.g. Refs. [G5L EH ETJ). 
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The functions F, in Eqs. (|2|)-(|6]) represent form factors for virtual mesons at the 
different vertices % (i — 7T, rj, f\) and for each vertex they are taken in the monopole form 

a 2 _ m i 

m) - ^, (8) 

where Aj is a cut-off parameter. For the 'effective' tt exchange we use the coupling 
constant fl NN /^ = 0.08 and cut-off parameter A nNN = 1.05 -j- 1.3 GeV. In the case of 
rj exchange we take g 2 NN /4n = 3, A r? 7 V v=l-5 GeV and <? aoT/7r =2.46 GeV, which results 
from the width T(a — > t]tt) = 80 MeV. 

The contribution of the /i exchange is calculated with g^NN = 11-2, Af lNN = 



1.5 GeV from Ref. |27[ and g ao / l7r =2.5. The latter value for g ao h^ corresponds to 



Ftot(fi) = 24 MeV and Br(fi — > a n) = 34%. The same parameters have been used in 



our previous study of a^ production in ttN — > a^N and pp — > da$ reactions |21[ . 

For the form factors at the ao/1 71 " ( a s well as cto^ 71 ") vertex factorized forms are 
applied following the assumption from Refs. [p8|, |2T 



Faofi n(ti,t 2 ) — F/iJVAT^i) F nNN (t 2 ), (9) 

where F^^^it), F n NN(t) are taken as in @. 

According to different versions of the Bonn potential the coupling constant 
daoNN/^ 71 can var y from 1.1075 to 2.67 p5[ 0. On the other hand, the unitary model 



for meson-nucleon scattering |30| gives a different range for this constant from 0.0026 to 
0.88. In the latter model the a only gives a contribution to the nr] background because 
there are no known resonances which decay to oqN. Since the model is extended only up 
to energies a/s < 1.9 GeV, which is below the a^ threshold, the meson-nucleon dynamics 
is not very sensitive to the oqNN coupling. We note that a small value of g^NN/^ 71 
certainly contradicts the experimental values of Br(pp —>■ a Tr) = 0.69 ± 0.12 |3l| and 
Br(pp — > oquj) = 0.354 ± 0.028 [^| , which are quite large (see e.g. Refs. []26, 



Having in mind these considerations we take (as well as in Ref. [^T]) the minimal value 
suggested by the Bonn potential g ao NN — 3.7. This value is not very different from the 



upper value of 3.33 given by the model of Ref. [|30| . 

Another problem is the treatment of a virtual nucleon. In this case - instead of 
the product of two monopole form factors (at the oqNN and ttNN vertices) - we use a 
dipole-like form factor, 

Ms) = A4 , ^ -TTo , (10) 



M + 



N -r \a - '"-at 



m 



which is normalized at s — m 2 and has the same asymptotics at large s (positive or 
negative) as Fi(s)Fj(s). 

There are a couple of arguments in favour of using the form factor ([T^) for virtual 



nucleons instead of those which are applied for virtual mesons. In the t-channel graph 
in elastic NN scattering the value of t is negative and the monopole form factor F n as 
given by Eq. (|8|) does not have a singularity in the physical region and decreases with 
t. For the s-channel graph with a nucleon exchange in the 7rN — > a N amplitude the 
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value of s is positive in the physical region and the conventional form factor 



A 2 



m 



N 



A%-s 

may have even a pole in the physical region (this happens for An = 2 GeV, which is used 
in the Bonn potential for a virtual do). This undesirable property is absent in the form 
factor (|TUD, where we consider the cut-off An as a free parameter. In our previous work 
f2~If we fixed An in the interval 1.2-1.3 GeV using experimental data on the differential 
cross section of the reaction pp — >■ da£ at pi a t, = 3.8 -f- 6.3 GeV/c JEJ; in this study we 
take A N = 1.24 GeV as an average value (see also the discussion in Section 4). 

We recall that the functional form of the nucleon form factor given by fllQl) was 
used in many papers, where meson production in irN, jN and NN collisions has been 
discussed (see e.g. |29], |3(J, |34|, |35], ^6[ [37| and references therein). 



The total cross section for a production in the isospin reaction j is given as 
the coherent sum of the amplitudes fll]) over all channels (a = s(N),u(N),t(fi),t(r])) 
integrated over phase space 

2 



a 



OQ » 



dE c dq dcos9 q d(p q 



2% 4 p aV /i 



J2Ml a) [ab;cd] 



111 



Here s = (p a + pb) 2 is the total energy of the NN system squared, E c and q are the 
energy of the outgoing nucleon and a$ meson, respectively. 6 q is the polar angle of the 
3-momentum of the ao-meson q in the cms of the initial nucleons defined as 9 q = q7p a) 
while ip q is the azimuthal angle of q in the cms. 



As shown in the analysis in Ref. [J2TJ the contribution of the ^-exchange to the 
amplitude irN — > a N is small. Note that in Ref. |38| only this mechanism was taken 
into account for the reaction pn — > ppa^ . Here we also include the r/-exchange because 
it might be noticeable in those isospin channels where a strong destructive interference 
of u- and s-channel terms can occur (see below). 



3. The Reggeon exchange model 



Here as in Ref. |2T| we also use the Regge-pole model for the amplitude ttN —>■ aoN 
as developed by Achasov and Shestakov [ 13] . The s-channel helicity amplitudes for the 
reaction n~p — » a^n in this approach can be written as 

B( 8 ,ty 



MvxA^P 



a Q n) 



u\>(p' 2 . 



-A(s,t) + (pi +Pi) Q 7o 



lsu x M%) 



where the invariant amplitudes A(s,t) and B(s,t) do not contain kinematical 
singularities. The relations between the invariant and s-wave helicity amplitudes are 
given by 

e 



M 



++ 



M 4 



-M_ 



cos • 



Ays, t)\/ —t m i n — B{s, t)yj —t max s 



cos ■ 



A(S, t)y/-t max ~ B(s, t)y/-t min S 



(13) 
(14) 
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where 9 is the cm. scattering angle, while t 
180°, respectively. 



, ,,,,„ and t max are the values of t at 0=0° and 



In the model of Ref. [15| the s-channel helicity amplitudes are expressed through 
the b\ and the conspiring p 2 Regge trajectories exchange as follows 

7 P2 (t)exp 



M. 



++ 



.71 . . 

-i-a P2 (t) 



(15) 



M 4 



\J (tnin - t) / S 7 bl (t) % exp 



,7T 



s > 



(16) 



As in Ref. [^T] we take the meson Regge trajectories in linear form atj(t) = <x,(0) + o^(0)£ 
with a bl {0) ~ -0.37, a P2 (0) ~ -0.6 and a^(0)= a' p2 {0) = 0.9 GeV~ 2 . The residues are 
parametrized in a convential way, 7 P2 (t) = 7 P2 (0) ex P(fyo2*)) 7&i(0 = 7&i(0) ex P(^6i^); 
all parameters were taken the same as in Ref. |21| . They correspond to two fits of the 
Brookhaven data on da/dt at 18 GeV/c ]39| found by Achasov and Shestakov fl5 ]: a) 
with pure p 2 contribution and b) with combined p 2 + &i contribution. 

The invariant amplitude corresponding to the diagram of Fig. Q can be written as 

B(a,t)' 



■Ml egge [ab;cd] 



u{p c 



■A(s, t) + (p ao +p d - p b ) a ^ a - 



75«(Pa) 



x «(Pd)75«(Pfc) x R(Pb,Pd)- 



(17) 



4. The reaction NN -> iVA^a r 



In order to demonstrate the sensitivity of the effective OPE model to the cut-off 
parameter A n NN used in the ttNN vertices we show in Fig. [3] the total cross section for 
the reaction pp — > pna^ for u(N) and t(fi) channels as a function of the excess energy 



Q = \/s - y/s~o~, where ^/s^ 



711 



«<> 



+ 2mN, calculated for different cut-off parameters. 



The dotted lines correspond to A wNN = 0.8 GeV, the solid lines show the result for 
A n NN = 1-05 GeV whereas the dashed lines indicate A^iv = 1-3 GeV. The results for 
A n NN = 1-3 GeV and 0.8 GeV differ by a factor of ~ 5. For our subsequent calculation 
we choose A^nn = 1-05 GeV while keeping the uncertainty on A-^nn in our 'effective' 
approach in mind. 

Since we have two nucleons in the final state it is necessary to take into account 
their final-state-interaction (FSI), which has some influence on meson production near 
threshold. For this purpose we adopt the FSI model from Ref. 0] based on the 
(realistic) Paris potential. We use, however, the enhancement factor Fnn{qnn) - as 
given by this model - only in the region of small relative momenta of the final nucleons 
Qnn < <7o? where it is larger than 1. Having in mind that this factor is rather uncertain 
at larger qNN, where for example contributions of nonnucleon intermediate states to the 
loop integral might be important, we assume that F NN (q NN ) = 1 for q NN > q Q . 

In Fig. | we show the FSI effect on the total cross section for the reactions pp — > ppa® 
(upper part) and pp — ► pna^ (lower part) for u(N), s(N) and t(/i) channels. The solid 
lines show the calculation without FSI whereas the dashed lines indicate the results with 
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FSI. As seen from Fig. [|, the FSI effect is stronger for pn than for pp in the final state 
due to the Coulomb repulsive interaction in the pp system and the isospin dependence 
of the NN interaction at small relative momenta. 

The results of our calculations for the total cross sections with FSI for the different 
isospin reactions are presented in Figs. ||, |6] as a function of Q = s/s — ■y/so- In Fig. ^| 
we show the total cross section for the pp reactions - pp — > ppa$ (upper part) and 
pp — > pna~Q (lower part), whereas in Fig. ^ we display the results for the pn reactions 
- pn —>■ ppaQ (upper part) and pn — > pnaQ (lower part). The solid lines with full dots 
and with open squares (r.h.s.) represent the results within the p 2 and (p2,bi) Regge 
exchange model. The short dotted lines (l.h.s.) corresponds to the t(j\) channel, the 
dotted lines to the t(n) channel, the dashed lines to the u(N) channel, the short dashed 
lines to the s(N) channel. The dashed line in the right upper part of Fig. |5| is the 
incoherent sum of the contributions from s(N) and u(N) channels (s + u). 

As seen from Figs. |5| and || the u- and s-channels give the dominant contribution; 
the t(/i) channel is small for all isospin reactions. For the reactions pp — > pna^, 
pn — > ppaQ and pn —>■ pna® the Regge exchange contribution (extended to low energies) 
becomes important and for the pn — > pna® reaction this contribution is even dominant 
near threshold. For the pp — > ppa$ channel the Regge model predicts no contribution 
from p 2 and P2,b\ exchanges due to isospin arguments (i.e. the vertex with a coupling 
of three neutral components of isovectors vanishes); thus only s-, u- and t(/i)- channels 
are plotted in the upper part of Fig. |[ 

Here we have to point out the influence of the interference between the s- and u- 
channels. According to the isospin coefficients from the OPE model presented in Table [I], 
the phase (of interference a) between the s- and u- channels M.1t N \ + exp(—ia)A4^ N ^ is 
equal to zero, i.e. the sign between Ai^r N \ and M*rm is 'plus'. The solid lines in Figs. ^, 
^ indicate the coherent sum of s(N) and u(N) channels including the interference of 
the amplitudes (s + u + int.). One can see that for pp — > pna^, pn — > ppa^ and 
pn — > pna$ reactions the interference is positive and increases the cross section, whereas 
for the pp —>■ ppa$ channel the interference is strongly destructive since we have identical 
particles in the initial and final states and the contributions of s- and -u-channels are 
very similar. 

Here we would like to comment about an extension of the OPE (one-pion- exchange) 
model to an OBE (one-boson-exchange) approximation, i.e. accounting for the exchange 
of a, p, u, ... mesons as well as for multi-meson exchanges. Generally speaking, the total 
cross section of Oq production should contain the sum of all the contributions: 

a(NN -> NNa ) = E^-, 

where j = K,a,p,u).... Depending on their cut-off parameters the heavier meson 
exchanges might give a comparable contribution to the total cross section for a 
production. An important point, however, is that near threshold (e.g. Q < 0.3 GeV 
) the energy behaviour of all those contributions is the same, i.e. it is proportional 
to the three-body phase space o~j ~ Q 2 (when the FSI is switched off and the narrow 
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resonance width limit is taken). In this respect we can consider the one-pion exchange 
as an effective one and normalize it to the experimental cross section by choosing an 
appropriate value of A w . The most appropriate choice for A w is about l-hl.3 GeV. 
Another question is related to the isospin of the effective exchange. As it is known from 
a serious of papers on the reactions NN — > NNX,X = n,r]',oo,(j) near threshold the 
most important contributions to the corresponding cross sections comes from n and p 
exchanges (see e.g. the review ]^lj and references therein). In line with those results 



we assume here that the dominant contribution to the cross section of the reaction 
NN — > NNa comes also from the isovector exchanges (like ir and p). In principle, it 
is also possible that some baryon resonances may contribute. However, as mentioned 
above, there is no information about resonances which couple to the OqN system. Our 
assumptions thus enable us to make exploratory estimates of the ao production cross 
section without introducing free parameters that would be out of control by existing 
data. The model can be extended accordingly when new data on the a production will 
be available. 

Another important question is related to the choice of the form factor for a 
virtual nucleon, that - in line with the Bonn-Julich potentials - we choose as given 
by (|10D, which corresponds to monopole form factors at the vertices. In the literature, 
furthermore, dipole-like form factors (at the vertices) are also often used (cf. Refs. 
|29| , |30| , [34], ^, |36|). However, there are no strict rules for the 'correct' power of the 
nucleon form factor. In physics terms, the actual choice of the power should not be 
relevant; we may have the same predictions for any reasonable choice of the power if the 
cut-off parameter A N is fixed accordingly. Note, that A^r may also depend on the type 
of mesons involved at the vertices. Therefore, we can not simply employ the parameters 
from Refs. [35| , J53 or others in case of the a problem. 



In our previous work [21] we have fixed A^v for the monopole related form factor (fLOf) 
in the interval 1.2-1.3 GeV fitting the forward differential cross section of the reaction 
pp — > da~Q from ||33| . On the other hand, the same data can be described rather well using 
a dipole form factor (at the vertices) with An =1.55-1.6 GeV (cf. Fig. [?]). If we employ 
this dipole form factor with A^ =1.55-1.6 GeV in the present case we obtain practically 
identical predictions for the cross sections of the channels pp — > pna^ , pn — > pna®, 
pn — > ppaQ , where the w-channel mechanism is dominant and u — s interference is not 
too important. In the case of the channel pp — > ppa® we obtain cross sections by up to a 
factor of 2 larger for the dipole-like form factor in comparison to the monopole one. This 
is related to the strong destructive interference of the s and u exchange mechanisms, 
which slightly depends on the type of form factor used. However, our central result, 
that the cross section for the pna$ final channel is about an order of magnitude higher 
than the ppa$ channel in pp collisions, is robust (within less than a factor of 2) with 
respect to different choices of the form factor. 

As seen from Figs. [], [], we get the largest cross section for the pp — > pna^ isospin 
channel. For this reaction the it-channel gives the dominant contribution, the s-channel 
cross section is small such that the interference is not so essential as for the pp — > ppa$ 
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reaction. 

The result within the Regge model is shown in Fig. ^ for the reactions pp —>■ pna$, 
pn — > pna^ (upper part) and pn — > ppa^ (lower part) in a wide energy regime 
for Q =1 MeV-i-10 GeV. The total cross section is calculated with the p 2 (dashed 
lines) and (p 2 ,bi) (solid lines) Regge trajectories (with FSI) for a cut-off parameter 
^■kNN = 1-05 GeV. In order to show the influence of the cut-off parameter A nNN in the 
Regge model we present in the lower part of Fig. || also the results for A nNN = 1.3 GeV 
(dotted line for p 2 exchange and the dot-dashed line for the (p 2 , &i) trajectory). Changing 
the cut-off A^n from 1.05 to 1.3 GeV gives a factor ~ 2 in the total cross section similar 
to the results within the effective Lagrangian model (cf. Fig. |||). 

As it was already discussed in our previous study [^l] an effective Lagrangian model 
cannot be extrapolated to high energies because it predicts the elementary amplitude 
it N — > aoN to rise fast. Therefore, such model can only be employed not far from the 
threshold; at larger energies it has to be unitarized. On the other hand, the Regge 
model is valid at large energies and we have to worry, how close to the threshold we can 
extrapolate corresponding amplitudes. According to duality arguments one can expect 
that the Regge amplitude can be applied at low energy, too, if the reaction ttN — > aoN 
does not contain essential s-channel resonance contributions. In this case the Regge 
model might give a realistic estimate of the ttN — > oqN amplitude even near threshold. 

Anyway, as we have shown in our previous paper [2Jj the Regge and w-channel model 



give quite similar results for the n~p — > a^n cross-section in the near threshold region; 
some differences in the cross sections of the reactions AW — > NNao - as predicted by 
those two models - can be attributed to differences in the isospin factors and effects of 
AW antisymmetrization which is important near threshold (the latter was ignored in 
the Regge model formulated for larger energies). 

5. The reaction AW -* NNa -> NNKK 

5.1. The KK and nn decay channels of the a (980) 

The ao(980) meson production has not yet been measured in AW — * NNao reactions. 
There are only a fewpp — > ppKK and pp — > pnKK experimental data points. Therefore, 
it is important to analyse a possible resonance contribution to KK production in the 
reactions AW — > NNX, using the calculated AW — > NNa amplitudes and the 
experimental fits obtained for the a^ resonance mass distribution in the KK decay 
channel. 

The amplitude for the ao(980) decays into KK and 7rrj modes can be parametrized 
by the well-known Flatte formula ||2| which satisfies both requirements of analyticity 
and unitarity for the two-channels ttt] and KK. 

In the case of the ao(980) resonance the mass distribution of the final KK system 
can be written as a product of the total cross section for a® production (with the 'running' 
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mass M) in the NN — > NNoq reaction ( |TT] ) and the Flatte mass distribution function 

lw {s,M) = a ao {s,M) C F {M2 _ Ml)2 + M2A{M) (18) 

with the total width T tot (M) = T aoKR (M) + T ao7Trt (M). The partial widths 

r ao ,r,(M) = Ql^^p (19) 

are proportional to the decay momenta in the center-of-mass (in case of scalar mesons), 

[(M 2 - (m K + m R ) 2 )(M 2 - (m K - m R ) 2 )] 1/2 



Qkk 

Qtttj ~ 



2M 
[(M 2 - (m* + m v ) 2 )(M 2 - (m n - m,) 2 )] 1/2 



2M 

for a meson of mass M decaying to KK and tttj, correspondingly. The branching ratios 
Br(a Q — > i^^) and Br(ao — > 7rr]) are given by the integrals of the Flatte distibution 
over the invariant mass squared dM 2 = 2MdM: 

7 dM 2 M C F M H r„„(M) 
*(<•»—>> = / (M* - Ml? + MITUM) (21) 

^r^ dM2MC F M R T ao7TV (M) 



(M 2 - M 2 R - M R T aoKR (M)) 2 + M 2 r 2 om? (M)- 

The parameters CF,g KR ,g nv have to be fixed under the constraint of the unitarity 
condition 

Br{a -> ATA") + 5r(a ->• 7177) = 1 . (22) 

Choosing the parameter r = r ao7rr) (M#) in the interval 50 -j- 100 MeV as given by 
the PDG [f43f , one can fix the coupling g nv according to ([T9|). In Ref. |45| a ratio of 
branching ratios has been reported, 

Brian —>■ KK) , . 

r(a 980 = — V -^ = 0.23 ± 0.05, 23 

Br{a — > 7T77J 



for m ao = 0.999 GeV, which gives Br(ao — > KA') = 0.187. In another recent study jP 
the WA102 collaboration reported the branching ratio 

T(a -> irir)/r(a -> Tir]) = 0.166 ± 0.01 ± 0.02, (24) 

which was determined from the measured branching ratio for the /i(1285)-meson. In 



our present analysis we use the results from f4q| , however, keeping in mind that this 
branching ratio Br(ao —>■ KK) more likely gives an 'upper limit' for the a — > KK 
decay 
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Thus, the two other parameters in the Flatte distribution Cp and g ao KR can 
be found by solving the system of integral equations, for example, Eq. ( p0|) for 
Br(ao — > KK) = 0.187 and the unitarity condition (|22|). For our calculations we 
choose either T ao7Tr) (Mn) = 70 MeV or 50 MeV, which gives two sets of independent 
parameters CF,g ao KKi9a TTri for a fixed branching ratio Br(a — > KK) = 0.187: 

set I (r aom? = 70 MeV) : (25) 

g aoKR = 2.297, g aonv = 2.189, C F = 0.365 
set 2 (T aowv = 50 MeV) : (26) 

g aoKR = 1.943, 9atm = 1.937, C F = 0.354. 

Note, that for the K + K~ or K°K° final state one has to take into account an isospin 
factor for the coupling constant, i.e. g ao K+K- — g ao K°K = 9a KR/V^, whereas 

g ao K+K° — g ao K-K° — g ao KK- 

5.2. Numerical results for the total cross section 

In the upper part of Fig. ^ we display the calculated total cross section (within parameter 
set 1) for the reaction pp —>■ pna^ — > pnK + K° in comparison to the experimental 
data for pp — > pnK + K° (solid dots) from Ref. |^|| as a function of the excess energy 
Q = y/s— y^so. The dot-dashed and solid lines in Fig. |9| correspond to the coherent sum 
of s(N) and u(N) channels with interference (s + u + int.), calculated with a monopole 
form of the form factor (|To|) with A^v = 1-24 GeV and with a dipole form of (0) with 
Ajv = 1-35 GeV, respectively. We mention that the latter (dipole) result is in better 
agreement with the constraints on the near-threshold production of a in the reactions 
7i~p — > K~K°p and tt + p —>■ K + K°p [fi7 |. In the middle part of Fig. |] the solid lines 



with full dots and with open squares present the results within the p 2 and (p 2 , &i) Regge 
exchange model. The short dashed line shows the 4-body phase space (with constant 
interaction amplitude), while the dashed line is the parametrization from Sibirtsev et 
al. f4"5f . We note, that the cross sections for parameter set 2 are similar to set 1 and 



larger by a factor ~ 1.5. 

In the lower part of Fig. ^| we show the calculated total cross section (within 
parameter set 1) for the reaction pp — > ppaf^ — > ppK + K~ as a function of Q = 
V$ ~ V^o ^ n comparison to the experimental data. The solid dots indicate the data 
for pp — > ppK°K° from Ref. [HB], the open square for pp — > ppK + K~ is from the 



DISTO collaboration [[4J| and the full down triangels show the data from COSY- 11 |5Q 

For the pp — > ppa$ — > ppK + K~ reaction (as for pp — > ppa$) there is no contribution 
from meson Regge trajectories; s- and w-channels give similar contributions such that 
their interference according to the effective OPE model (line s + u + int.) is strongly 
destructive (cf. upper part of Fig. ||). The t(f\) contribution (short dotted line) is 
practically negligible, while the t(^)-channel (dotted line) becomes important closer to 
the threshold. 

Thus our model gives quite small cross sections for a° production in the pp — > 
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ppK + K~~ reaction which complicates its experimental observation for this isospin 
channel. The situation looks more promising for the pp — > pna^ — > pnK + K° reaction 
since the a^ production cross section is by an order of magnitude larger than the a$ 
one. Moreover, as has been pointed out with respect to Fig. |], the influence of the 
interference is not so strong as for the pp — > ppaf^ — > ppK + K~ reaction. 

Here we stress again the limited applicability of the effective Lagrangian model 
(ELM) at high energies. As seen from the upper part of Fig. [|, the ELM calculations 
at high energies go through the experimental data, which is not realistic since also other 
channels contribute to K + K° production in pp reactions (cf. dashed line from Ref. ||48| ). 
Moreover, the ELM calculations are higher than the Regge model predictions which 
indicates, that the ELM amplitudes at high energies have to be reggeized or unitarized. 

5. 3. Numerical results for the invariant mass distribution 

As follows from the lower part of Fig. [], the oq contribution to the K + K~ production 
in the pp — ► ppK + K~ reaction near the threshold is hardly seen. With increasing 
energy the cross section grows up, however, even at Q — 0.111 GeV the full cross 
section with interference (s + u + int.) gives only a few percent contribution to the 
0.11 ± 0.009 ± 0.046 fjb 'nonresonant' cross section (without <fi — > K + K~) from the 
DISTO collaboration |9|. 

To clarify the situation with the relative contribution of a° to the total K + K~ 
production in pp reactions we calculate the K + K~ invariant mass distribution for the 
pp — » ppK + K~ reaction at p^ = 3.67 GeV/c, which corresponds to the kinematical 
conditions for the DISTO experiment ||49|| . The differential results are presented in 
Fig. [10]. The upper part shows the calculation within parameter set 1, whereas the lower 
part corresponds to set 2. The dot-dashed lines (lowest curves) indicate the coherent 
sum of s(N) and u(N) channels with interference (s + u + int.) for the a^ contribution. 
However, one has to consider also the contribution from the /o scalar meson, i.e. the 
pp — > ppfo — > ppK + K~ reaction. The /o production in pp reactions has been studied 
in detail in Ref. plf . Here we use the result from this previous work [|51j] and show 



in Fig. [H] the contribution from the fo meson calculated with parameter set A from 
Ref. [^TJ as the solid line with open circles (fo). 

We find that when adding the /o contribution to the phase-space of nonresonant 
K + K~ production (the short dotted lines in Fig. [ll]) and the contribution from decays 
(resonance peak around 1.02 GeV), the sum (solid) lines almost perfectly describe the 
DISTO data. This means that there is no visible signal for an a[] contribution in the 
DISTO data according to our calculations while the /o meson gives some contribution 
to the K + K~ invariant mass distribution at low invariant masses M, that is ~ 12% 



of the total 'nonresonant' cross section from the DISTO collaboration |49|]. Thus the 
reaction pp — > pnK + K° is more promising for a^ measurements as it has been pointed 
out in the previous subsection. 

For an experimental determination of the a$ we present the invariant mass 
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distribution of K + K° in the reaction pp — > pnK + K° at different Q (solid lines) in 



Fig. [[J]. The dashed lines show the invariant mass distributions for 'background' (i.e. 
according to phase space with constant interaction amplitude) under the assumption 
that the integrals below the solid and dashed lines are the same for each Q. We see 
that the shape of the solid and dashed lines are practically the same for Q < 50 MeV. 
Noticeable differences between the lines can be found for Q > 100 MeV. This means that 
a separation of the resonance contribution from the background very close to threshold 
can be done only in the case when the background is small or very well known. 

6. Conclusions 

In this work we have estimated the cross sections of a^ production in the reactions 
pp — > ppa^, pp — >■ pnaQ, pn — > ppa^ and pn — > pnag near threshold and at medium 
energies. Using an effective Lagrangian approach with one-pion exchange we have 
analyzed different contributions to the cross section corresponding to t-channel diagrams 
with ^(550)- and /i(1285)-meson exchanges as well as s and u-channel graphs with an 
intermediate nucleon. We use the same parameters as in our previous paper where we 
describe rather well the Berkeley data [|33| on the reaction pp — > da^. 



We additionally have considered the t-channel Reggeon exchange mechanism with 
parameters normalized to the Brookhaven data for ir~p — > a$p at 18 GeV/c |39|| . These 
results have been used to calculate the contribution of ao mesons to the cross sections 
of the reactions pp — ► pnK + K° and pp — > ppK + K~ . Due to unfavourable isospin 
Clebsh-Gordan coefficients as well as rather strong destructive interference of the s- and 
w-channel contributions our model gives quite small cross sections for a® production 
in the pp — > ppK + K~ reaction. However, the a$ production cross section in the 
pp — > pnaQ — > pnK + K° reaction should be larger by about an order of magnitude. 
Therefore the experimental observation of a$ in the reaction pp — > pnK + K° is much 
more promising than the observation of a® in the reaction pp — > ppK + K~ . We note in 
passing that the tttj decay channel is experimentally more challenging since, due to the 



larger nonresonant background | 52| , the identification of the r/-meson (via its decay into 
photons) in a neutral-particle detector is required. 

We have also analyzed invariant mass distributions of the KK system in the reaction 
pp — > pNao — > pNKK at different excess energies Q not far from threshold. Our analysis 
of the DISTO data on the reaction pp — > ppK + K~ at 3.67 GeV/c has shown that the 
aQ-meson is practically not seen in da/dM at low invariant masses, however, the /o- 
meson gives some visible contribution. In this respect the possibility to measure the 
aj" meson in da/dM for the reaction pp — > pnK + K° (or — > dK + K°) looks much more 
promising not only due to a much larger contribution for the aj", but also due to the 
absence of the /o meson in this channel. 

Experimental data on ao production in NN collisions are practically absent (except 
of the ao observation in the reaction pp — > dX [33|). Such measurements might give new 



information on the ao structure. According to Atkinson et al. |53| a relatively strong 
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production of the a$ (the same as for the &i(1235)) in non-diffractive reactions can be 
considered as evidence for a qq state rather than a qqqq state. For example the cross 
section of a® production in jp reactions at 25-50 GeV is about 1/6 of the cross sections 
for p and u production. Similar ratios are found in the two-body reaction pp — > dX at 
3.8-6.3 GeV/c where a{pp — ► da$) = (1/4 -f- l/6)a(pp — > dp + ). 

In our case we can compare a and uj production. Our model predicts a{pp — > 
pna^) = 30 -f- 70/xb at Q ~ 1 GeV (see Fig. ||) which can be compared with 
aijpp — >■ ppu) ~ 100 — 200/ib at the same Q. If such a large cross section could be 
detected this would be a serious argument in favour of the qq model for the a . 

To distinguish between the threshold cusp scenario and a resonance model one 
can exploit different analytical properties of the a production amplitudes in those 
approaches. In case of a genuine resonance the amplitude of r)7r and KK production 
through the a has a pole and satisfies the factorization property. This implies that 
the shapes of the invariant mass distributions in the rjir and KK channels should not 
depend on the specific reaction in which the ao resonance is produced (for Q > Y tot ). On 
the other hand, for the threshold cusp scenario the a$ bump is produced through the irrj 
final state interaction. The corresponding amplitude has a square root singularity and 
in general can not be factorized (see e.g. Ref. |4(| were the factorization property was 



disproven for pp FSI in the reaction pp — > ppM). This implies that for a threshold bump 
the invariant mass distributions in the nir and KK channels are expected to be different 
for different reactions and will even depend on kinematical conditions (i.e. initial energy 
and momentum transfer) at the same exess energy, e.g. Q ~ 1 GeV. 
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Reaction j (mechanism a) €j( a )[ab; cd] £]( a )[ab; dc] £]( a )[ba;dc] £]( a )[ba; cd] 
PP^Ppa° Q {t{ri),t{h)) +1/V2 -1/V2 +T772 -1/V2 

(s(N)) +1/V2 -1/V2 +1/V2 -1/V2 

(u(N)) +1/V2 -1/V2 +l/v^ -1/V2 

Regge 



pp^pna^ (t(rj),t{fi)) 


-V2 








+V2 


(s(N)) 





+v/2 


-v/2 





(u(N)) 


+2v/2 


-v/2 


+v / 2 


-2^2 


Regge 


-1 


+ 1 


-1 


+1 


pn ^ppa^" (t(r/),t(/i)) 


+1 


-1 








( S (iV)) 


-2 


+2 


-1 


+1 


(u(JV)) 








+ 1 


-1 


Regge 


+1/V2 


-1/V2 


-1/V2 


+1/V2 


pn ^pnag (t(rj),t(fi)) 


-1 





+ 1 





(s(N)) 


-1 


-2 


+1 


+2 


(u(N)) 


-1 


+2 


+1 


-2 


Regge 





+\/2 





-72 



Table 1. Coefficients in Eq. rtjl) for different mechanisms of the pp — > ppap, 
pp — ► ppaj, pn — ► ppaQ and pn — * pna[J reactions. 
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Figure 1. Diagrams for ao production in the reaction NN — > ciqNN near threshold 
as considered in the present study. 
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Figure 2. The diagram for ao production in the reaction NN — > NNao within the 
Rcgge exchange model. 
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Figure 3. The total cross section for the reaction pp — ► pna^ for u(N) and i(/i) 
channels as a function of the excess energy Q = i/s — ^/sq for different cut-off 
parameters A^nn = 0.8 GeV (dotted lines), A^nn — 1-05 GeV (solid lines) and 
A^nn = 1.3 GeV (dashed lines). 
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Figure 4. The total cross section for the reactions pp — > ppa® (upper part) and 
pp — > pna,} (lower part) as a function of the excess energy Q = y's — -^so f° r u(N), 
s(N) and t(/i) channels calculated without FSI (solid lines) and with FSI (dashed 
lines). 
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Figure 5. The total cross sections for the reactions pp — > ppa^ (upper part) and 
pp — ► pna,Q (lower part) as a function of the excess energy Q = y/s — y/so calculated 
with FSI. The short dotted lines (l.h.s.) corresponds to the £(/i) channel, the dotted 
lines to the t{rj) channel, the dashed lines to the u(N) channel, the short dashed lines 
to the s(N) channel. The dashed line (upper part, r.h.s.) is the incoherent sum of 
the contributions from s(N) and u(N) channels (s + u). The solid lines indicate the 
coherent sum of s(N) and u(N) channels with interference (s + u + int.). The solid 
lines with full dots and with open squares (lower part, r.h.s.) present the results within 
the P2 and (p2, &i) R-egge exchange model. 
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Figure 6. The total cross sections for the reactions pn — > ppa^ (upper part) and 
pn — » pna[J (lower part) as a function of Q = \/s — y/so calculated with FSI. The 
assignment of the individual lines is the same as in Fig. ||. The results from the 
effective OPE model are shown on the l.h.s. while those from the pi and p^b\ Regge 
approach are displayed on the r.h.s. 



Production of clq -mesons in pp and pn reactions 



22 



a 



10 



3 pp^da^ 



10 



10 




10 



mon, A N =1.2 GeV 
dip,A N =1.2GeV 
dip, A N =1.55GeV 



^=1.05 GeV 



mon, A N =1 .24 GeV 

- - dip,A N =1.24GeV 
dip, A N =1.6GeV 



10 



10 



10 



p lab -3.29 [GeV/c] 



Figure 7. Forward differential cross section of the reaction pp — » da^ as a function of 
(piab — 3.29) GeV/c. The bold and thin solid curves are calculated at A T 7VA r =l-05 and 
1.3 GeV, respectively. The solid curves correspond to a monopole nucleon form factor 
with An= 1.2 (thin) and 1.24 GeV (bold). The long-dashed and short-dashed curves 
are calculated using the dipole nucleon form factor for different values of An as shown 
in the figure. The experimental data are taken from Ref. J33J. 
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Figure 8. The total cross sections for the reactions pp — > pna^ , pn — > pna^ (upper 
part) and pn — ► ppa^ (lower part) as a function of Q — \fs — ^/sq calculated within the 
pi (dashed lines) and (p2,b\) (solid lines) Regge exchange model (with FSI) for cut-off 
parameters A^nn = 1-05 GeV. The dotted and dot-dashed lines in the lower part show 
the results for A^nn = 1.3 GeV within the pi and (p2,bi) exchanges, respectively. 
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Figure 9. Upper part: the calculated total cross section (within parameter set 1) 
for the reaction pp — > pna^ — » pnK + Ko in comparison to the experimental data for 
pp — > pnK + K (solid dots) from Ref. Q as a function of Q — y/s — ^/sq. The dot- 
dashed and solid lines correspond to the coherent sum of s(N) and u(N) channels 
with interference (s + u + int.) calculated with a monopole form of the form factor 
@ with A N = 1.24 GeV and with a dipolc form of @ with A N = 1.35 GeV, 
respectively. Middle part: the solid lines with full dots and with open squares represent 
the results within the p2 and {p2,bi) Regge exchange model. The short dashed line 
shows the 4-body phase space (with constant interaction amplitude); the dashed line 
is the parametrization from Sibirtsev et al. |4q] . Lower part: the calculated total cross 
section (within parameter set 1) for the reaction pp — > ppa^ — > ppK + K~ as a function 
of Q = y/s — ,/sq in comparison to the experimental data. The solid dots indicate 
the data for pp — ► ppKoKo from Ref. pq |, the open square for pp — > ppK + K~ from 
Rcf. Est] ; the full down triangels show the data from Ref. pCf] . 
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Figure 10. The K + K~ invariant mass distribution for the pp — > ppK + K~ reaction 
at pi a b = 3.67 GeV/c. The short dotted lines indicate the 4-body phase space with 
constant interaction amplitude, the dot-dashed lines show the coherent sum of s(N) 
and u(N) channels with interference (s + u + int.). The solid lines with open circles 
correspond to the /o contribution from Ref. |£l[. The thick solid lines show the sum 
of all contributions including the decay <p ~ * K + K~ . The experimental data are taken 
from Ref. II. 
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Figure 11. The K + K° invariant mass distribution for the pp — * pnK + K° reaction 
at different Q = y/s — y/s . The solid lines describe the aj resonance contributions. 
The dashed lines show the invariant mass distributions for 'background' under the 
assumption that the integrals below the solid and dashed lines are the same for each 
Q. 



